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Fast Scatter Artifacts Correction for Cone-Beam CT 
without System Modification and Repeat Scan 

Wei Zhao, Jun Zhu, Luyao Wang 


Abstract —We provide a fast and accurate scatter artifacts 
correction algorithm for cone beam CT (CBCT) imaging. The 
method starts with an estimation of coarse scatter profile for 
a set of CBCT images. A total-variation denoising algorithm 
designed specifically for Poisson signal is then applied to derive 
the final scatter distribution. Qualitatively and quantitatively 
evaluations using Monte Carlo (MC) simulations, experimental 
CBCT phantom data, and in vivo human data acquired for a 
clinical image guided radiation therapy were performed. Results 
show that the proposed algorithm can significantly reduce scatter 
artifacts and recover the correct HU within either projection 
domain or image domain. Further test shows the method is robust 
with respect to segmentation procedure. 

Index Terms —Cone beam computed tomography (CBCT), 
scatter correction, denoising, polychromatic reprojection, spec¬ 
trum estimation. 


1. Introduction 

D ue to the increased axial coverage of multi-slice com¬ 
puted tomography (CT) and the introduction of flat- 
panel detector, the size of X-ray illumination flelds has grown 
dramatically, causing an increase in scatter radiation. For most 
of cases, artifacts are present when the model used in the 
reconstruction algorithm is not consistent with the projection 
data acquisition model. Existing reconstruction algorithms 
don’t model the scatter radiation, so scatter artifacts appear 
in the CT images. Typical scatter artifacts show as shading, 
streaks between high contrast objects, reduced contrast reso¬ 
lution and inaccurate Hounsfleld Units (HU) numbers. Scatter 
correction has been extensively studied in the past serval 
decades but a clinically sensible solution remains illusive. 
Current scatter correction methods can be roughly classifled 
into flve approaches: physical scatter rejection, analytical 
modeling, Monte Carlo (MC) simulation, primary modulation 
and scatter measurements |[T]. 

To develop a scatter artifacts correction method that is 
fast, accurate, and, ideally, has no need for extra scans, extra 
measurements or extra hardware support or modiflcations is 
the goal of this study. We provide a fast and accurate scatter ar¬ 
tifacts correction method for cone beam CT (CBCT) imaging. 
In the proposed method, the difference between the measured 
raw projection data and a polychromatic reprojection of a 
segmented image volume is used to generate a coarse estimate 
of scatter profile. Based on the coarse scatter estimation, a 
denoising algorithm specifically designed for Poisson signal 
which employs a total variation regularization term was then 
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applied to yield the final scatter estimation. By simply sub¬ 
tracting the final scatter estimate from the raw projection data, 
we can obtain the corrected projection data. The corrected 
image was then reconstructed using the corrected projection 
data. MC simulations, experimental phantom data and in vivo 
human data were used to validate and evaluate the proposed 
scatter correction method. 

II. Material and methods 

The measured raw projection data can be modeled as the 
summation of the primary and scatter signals. A general 
scatter artifact correction strategy is to estimate the scattered 
radiation distribution and then subtract it from the measured 
raw projection data, i.e. 

Ip{a,x) = I{a,x) - Is{a,x), (1) 

where Ip is the corrected projection data, namely, the estimated 
primary projection data, / is the raw data and Is is the scatter 
signal. Indices a and x stand for projection view angle and 
detector channel number respectively. For simplicity, we will 
drop the index x in the following discussions. 


A. Coarse scatter estimation 

The proposed method, as shown in Fig.[^ starts with an im¬ 
age segmentation of the scatter artifacts contaminated CBCT 
image volume, followed by a polychromatic reprojection of 
the segmented image volume using a given x-ray spectrum 
and known data acquisition geometry. The reprojected data is 
then subtracted from the measured raw projection data at each 
given view angle to generate a coarse estimate of the necessary 
scatter profile used for denoising algorithm. 

Based on the segmented image volume, the primary signal 
can be modeled as follows: 


Ip = N 
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where N is the total number of photons, Q{E) is the polychro¬ 
matic X-ray spectrum, r]{E) is the energy-dependent efficiency 
of the detector, which is simply considered as proportional 
to photon energy E because most CT scanners use energy- 
integrating detectors. Emax is the maximum photon energy 
of the spectrum. ji^E^s) is the energy-dependent linear at¬ 
tenuation coefficient and I is the propagation path length for 
each ray, and can be calculated using a GPU-based ray-tracing 
algorithm 0 0 By using these notations, the flood field /q 
can be modeled as follows: 

P E-max 
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Fig. 1. Flowchart of the proposed scatter correction method. The method can 
starts with either the raw projection data or an scatter contaminated volume 
as input. After an initial CT image reconstruction, the image is segmented 
into different components, based on which a polychromatic reprojection is 
performed with predetermined spectrum. The reprojection data is subtracted 
from the raw projection data to yield a coarse scatter estimate. The coarse 
scatter is then applied to a denoising algorithm which was specifically 
designed for Poisson signal to yield the final scatter estimate. The corrected 
projection data is obtained by subtracting the final scatter from the raw 
projection data and the corrected CT image can be generated using a standard 
filtered backprojection reconstruction algorithm. 


yield a smooth scatter distribution |[^. The denoising algorithm 
which is based on the total-variation regularization is aimed 
to solve the following optimization problem, 

/s = argminj^ J dx{Is - ijoglg) + ^ j dx\VIs\^. (6) 

The first term of I© is a data-fidelity term that is designed 
specifically for Poisson statistical signal and keep Ig close 
to the data Ig, while the second term is a total-variation 
regularization term to keep the solution Ig smooth, namely, 
being dominated with low frequency content. is a constant 
to determine the relative importance of the two terms. This 
objective function is convex and can be solved using an 
variational approach j^. 

It has to note that the proposed method can also be imple¬ 
mented in image-domain. For those cases when we only have 
images or volume, we need to estimate scatter artifact error 
images which can be added directly to the uncorrected images. 
To achieve this goal, a scatter projection error which 
can be added linearly in the logarithmic raw-data domain can 
be calculated. Because the tomographic reconstruction is a 
linear process and the order of summation and backprojection 
operations are interchangeable. 


After replacing N with Jq and substituting it into Q, the result 
describes the estimated primary signal with flood field Jq: 
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Subtracting the above estimated primary signal from the 
measured total projection data /, the coarse scatter signal Ig 
can be estimated as follows: 


Is = I- 
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Note that the segmented image volume is generated from 
a scatter contaminated reconstruction and the polychromatic 
x-ray spectrum used in the reprojection can be recovered 
using an indirect transmission measurement-based spectrum 
estimation method Q. 


B. Denoising the coarse scatter 

Since the coarse scatter estimate Ig is highly dependent on 
the segmentation procedure, it may yield inaccurate results, 
especially for low contrast objects that have similar attenuation 
properties as the background and for the edges of two neigh¬ 
boring materials. To compensate for the inaccuracy caused by 
the inaccurate segmentation, instead of regularizing the coarse 
scatter using a convolution-based scatter model in this 
study, we directly denoise the coarse scatter estimate using a 
statistical-based denoising algorithm. 

As we know, scatter radiation distribution is typically 
predominantly of low frequency content not only in spatial 
domain but also in temporal domain. A denoising technique 
specifically designed for Poisson signal was applied to Ig to 


C. Monte Carlo simulations 

To validate the proposed algorithm, an anthropomorphic^ 
thorax phantom and an abdomen phantom were used to 
generate MC simulation data with a Geant4-based MC simu¬ 
lation package GATE Q. To quantitatively and qualitatively 
investigate the effect of scatter artifacts reduction, both of 
the phantoms include a water insert at the central area and 
the water insert has three small low contrast inserts (adipose, 
breast and liver) and a bone insert. 

Primary projection data, scatter only projection data, and 
primary plus scatter projection data were extracted indepen¬ 
dently. For a faster simulation, a parallel geometry and a plane 
X-ray source (2D 320 x 120 mm^ rectangle source) were 
used in the MC simulation. The distance from the source 
to the center of rotation is 750 mm and the distance from 
the detector to center of rotation is 450 mm. A circular scan 
was simulated and a total of 360 projections per rotation are 
acquired over an angular range of 360°. The polychromatic 
X-ray source spectrum is 125 kVp and it is generated using 
Spektr software with 5 mm aluminum filtration. For each 
of the simulations, a total number of 3 x 10^^ photons were 
emitted. 

In a first experiment, we investigate scatter correction for 
the anthropomorphic thorax phantom and abdomen phantom 
in both projection domain and image domain. The primary 
image was used as reference images and the total images were 
applied to the proposed scatter correction algorithm. 

To test the robustness of method with respect to the segmen¬ 
tation, two components (lung and tissue), three components 
(lung, tissue and bone) and four components (lung, water, 
tissue and bone) segmentation were performed based on the 
uncorrected thorax phantom images respectively. 

^http://www.qrm.de/content/pdf/QRM-Thorax.pdf. 
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True Scatter Coarse Scatter Denoised Scatter 



Fig. 2. Sinograms of true scatter, coarse scatter, and denoised scatter for 
an anthropomorphic thorax phantom and an abdomen phantom. The coarse 
scatter which is the difference of the raw projection data and the polychromatic 
reprojection data have global scatter profiles but also contain errors which 
were partially generated from inaccurate segmentation. The denoised coarse 
scatter smooth these errors and fit the true scatter profiles well. 


D. Physical phantom experiments and patient study 

The proposed algorithm was also evaluated using experi¬ 
mental data of a CatphanSOO phantom (The Phantom Labo¬ 
ratory, Salem, NY) scanned using a CBCT on-board imager 
(Varian 2100EX System, Varian Medical Systems, Palo Alto, 
CA). A total of 678 projections were evenly acquired in a 360 
degree rotation with 2x2 binning. Both wide collimation and 
narrow collimation modes were applied with the same scan 
parameters. 

A set of in vivo pelvis data was also used to evaluate the 
proposed method. This data was acquired in half-fan mode 
using State-of-the-art of the kV CBCT imaging system of 
the Varian TrueBeam system(Varian Medical Systems, Palo 
Alto, USA). The tube potential was set to 125 kVp and 
678 projections were acquired on a flat-panel detector (Varian 
4030CB imager) with 1024 columns and 768 rows. 


Thorax Phantom Water Insert Abdomen Phantom Water Insert 



Fig. 3. Results of the scatter correction for the MC simulation data of the 
thorax phantom and the abdomen phantom within both projection domain 
and image domain. (First row) Results for the primary signal and they can be 
regarded as ground truth. (Second row) Results for the primary plus scatter 
signal (total signal) and scatter induced shading artifacts and streaks are clearly 
seen. (Third row) Results for the scatter corrected images within projection 
data domain. (Fourth row) Results for the scatter corrected images within 
image domain. Shading artifacts and streaks are significantly reduced after 
correction within both projection domain and image domain. Display window: 
[-1200HU, 500HU] for thorax phantom images and [-300HU, 300HU] for the 
water insert images and the abdomen phantom. 
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III. RESULTS 

A. Monte Carlo simulations 

Eig. 1^ shows the sinograms of the MC reference scatter, 
coarse scatter and the denoised scatter for the anthropomorphic 
thorax phantom and the abdomen phantom. Eig. shows the 
results of the scatter correction using the proposed method 
for the MC simulation data of the thorax phantom and the 
abdomen phantom within both projection domain and image 
domain. The primary images were reconstructed using primary 
projections and they were regarded as scatter-free images. 
Eig. 1^ depicts the HU numbers of breast, adipose, liver and 
bone inserts (shown in Eig.[^ of primary images, total images 
and scatter corrected images. 

The influence of different segmentation methods on the 
accuracy of the scatter correction for the thorax phantom is 
depicted in Eig. 


Fig. 4. Results of reconstruction values in HU numbers at different ROIs 
of the water insert of (a) the thorax phantom and (b) the abdomen phantom 
before and after correction within both projection domain and image domain. 
The HU numbers are successfully recovered after scatter correction using the 
proposed method. PD and ID stand for projection domain and image domain, 
respectively. 

B. Experimental phantom and patient study 

Eig. shows the CT ima^s and the corresponding differ¬ 
ence images of the Catphan®500 phantom with and without 
scatter correction. Eig. [7] shows CT images of an in vivo 
pelvis scan with the kilovoltage CBCT imaging system of 
the Varian TrueBeam system (Varian Medical Systems, Palo 
Alto, USA) without and with scatter correction. After scatter 
correction, the HU number of the dark region was successfully 
recovered and the missing anatomical structures are recovered. 
It took about 1 min to correct a typical Varian clinical dataset 
(512 X 512 X 81) using a NVIDIA GeEorce GTX 480 card. 
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Corrected 


Fig. 5. The influence of the segmentation methods on the accuracy of the 
scatter correction is shown for the thorax phantom with water insert in both 
projection domain (PD) and image domain (ID). (First row) Results for the 
four components segmentation with lung, water, tissue and bone. (Second 
row) Results for the three components segmentation with lung, tissue and 
bone. (Third row) Results for the two components segmentation with only 
lung and tissue. Display window: [-1200HU, 500HU] for the thorax phantom 
images and [-300HU, 300HU] for the water insert images. 


Narrow Image Wide Image Corrected (PD) Corrected (ID) 



Fig. 6. Catphan®500 phantom with and without scatter correction. The 
narrow collimation image is considered to be the scatter-free image. Shading 
and reduced HU are shown in the wide collimation image.After correction, 
scatter artifacts are greatly reduced and HU are successfully recovered for 
both projection domain (PD) and image domain methods (ID). The difference 
images show each image subtracted with the narrow collimation image. Note 
that the wide collimation scan and the narrow collimation scan are two 
independent scans and the registration can not be perfect. Display window: 
[-150 HU, 150 HU] for both the CBCT images and the difference images. 

IV. Discussion and conclusions 

It has to note that the selection of the denoising parameters 
depends on segmented image. Since the noise level of the 
MC simulation data is very high, a relatively large K and 
p (compared to the patient data) were used to smooth the 
segmentation error and to yield a denoised scatter that fit the 
true scatter well enough. To our belief, most of the clinical 
images have much lower noise level and the segmented images 
would be superior to the segmented images of the MC data in 
this study. Nevertheless, a larger K and (3 are recommended 
for badly segmented images to reinforce regularization. 

In summary, this work presents a novel strategy to estimate 
a coarse scattered radiation profile, based on which a denoising 
algorithm was applied to yield the final scatter estimation, we 
demonstrated that a significant increase in image uniformity 
and HU accuracy were achieved after correction. The correc¬ 
tion algorithm requires minimum computational cost, requires 


Fig. 7. Scatter correction for a pelvis CBCT scan using the kV imaging system 
of the Varian TrueBeam system. Axial view, coronal view and sagittal view 
are depicted respectively. The uncorrected images acquired on the commercial 
system show shading artifacts and part of the anatomic structure was missing. 
After scatter correction, the shading artifacts were successfully reduced and 
the missing anatomic structure can be clearly seen. Display window: [-200 
HU, 100 HU] for all images. 

no modifications to existing system hardware, and can also 
handle cross-scatter for dual source dual energy CT scanner. 
Furthermore, the method can also be applied to spiral cone- 
beam CT where the anti-scatter grid can be uninstalled to 
further reduce radiation dose. 
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